Disorder-induced noncollinear ferromagnetism in models for (III,Mn)V 

semiconductors 
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We study the ground state properties of kinetic-exchange models for (III,Mn) V semiconductors with 
randomly distributed Mn ions. Our method is embedded in a path integral spin-wave type formalism 
leading to an effective action for Mn spins only with full Matsubara frequency dependence. The zero- 
frequency contribution to this action is equivalent to static perturbation theory and characterizes the 
stability of a given spin configuration, while the component linear in frequency can be interpreted 
as the joint Berry phase of the Mn and carrier system. 

For simple parabolic-band carriers the collinear ferromagnetic state with all Mn spins in parallel is 
always stationary but generically unstable. This instability can be characterized in terms of inverse 
participation ratios and is due to long-ranged nonlocal spin fluctuations. We also present results for 
the ground state magnetization as a function of an external field. 

For carrier dispersions involving anisotropy induced by spin-orbit coupling the collinear state is not 
even stationary and therefore also not the ground state. This interplay between the anisotropy in 
the carrier system and the disorder in the Mn positions reflects recent findings by Zarand and Janko 
(Phys. Rev. Lett. 89, 047201 (2002)) obtained within the RKKY approximation. The stationarity 
of the collinear state (with the magnetization pointing in one of the cubic symmetry directions) is 
restored in the continuum or virtual crystal approximation where disorder is neglected. 



I. INTRODUCTION 

Since several years diluted ferromagnetic semiconduc- 
tors are in the focus of research in solid state physics; for 
reviews see 1-8 . This is on the one hand because these sys- 
tems provide interesting and challenging physical prob- 
lems on their own right. On the other hand this big deal 
of interest is due to possible applications of these ma- 
terials in the emerging field of spintronics 9 ' 10 since they 
offer the perspective of combining ferromagnetism with 
the readily tunable transport properties of semiconduc- 
tors. 

An important achievement in this field was the fabri- 
cation of diluted Mn-doped GaAs via low temperature 
molecular beam epitaxy by Ohno and collaborators in 
1996 11 . This material showed a Curie temperature of 
110K, a result which has been reproduced in the mean- 
time by several other groups. Moreover, very recently re- 
ports on Mn-doped (III,V) semiconductors having Curie 
temperatures of room temperature or higher have ap- 
peared, and also other combinations of magnetic ions 
and/or host materials look prospective 12-19 . 

These recent developments have also already generated 
a large amount of theoretical research on diluted ferro- 
magnetic semiconductors 5-8 ' 20-53 . Such studies include 
electronic structure calculations based on density func- 
tional techniques 8 ' 51-53 , and investigations using specific 
models for such systems 5-7 ' 20-50 . In particular, very re- 
cently a series of studies has appeared on disorder-related 
phenomena in different models for diluted ferromagnetic 
semiconductors 29-32 ' 36-48 . In the present work we em- 
ploy kinetic-exchange models for carrier-mediated ferro- 



magnetism occurring in Mn-doped (III,V) semiconduc- 
tors which have been the basis of a large body of previous 
work 5,20-34 . The general Hamiltonian reads 



TL = Hi 



E 



rf 3 rJ(f- Ri)s(f) ■ S7 



(1) 



These models show several features in accordance with 
experiments. They consist of a kinetic term Tlkm for free 
valence-band carriers (holes) whose spin density s(f) is 
antifcrromagnetically exchange-coupled to localized Mn 
spins Si of length S = 5/2 at locations Ri by an spatially 
extended coupling J(r ), which we take to be of the form 



J(f) 



J, 



pel 
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The regularization parameter ao reflects the spatial range 
of the exchange coupling 46 . In a minimal description of 
carrier-induced ferromagnetism the kinetic term repre- 
sents just a simple parabolic band characterized by an 
effective mass m*. Moreover we will also consider the 
case of a more realistic k ■ p Hamiltonian 54 ' 55 ' 26 describ- 
ing the valence band structure of (III,V) semiconductors 
such as GaAs. 

Differently from earlier work 20-27,31 we will not make 
use of the virtual crystal approximation, i.e. we will not 
approximate the localized Mn moments by a continuum 
but retain them as individual and randomly distributed 
spins. This additional feature of the model grossly en- 
riches its physical properties and gives rise to the oc- 
currence of noncollinear ferromagnetism 30 . This non- 
collincarity in the orientations of localized magnetic mo- 
ments even in the ground state of the system is very 
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likely to be an ingredient to the interpretation of recent 
experiments by Potashnik et aL 56 who found a strong de- 
pendence of the magnetic properties of Mn-doped GaAs 
on the annealing history of the sample. Similarly strong 
dependencies are also found in transport measurements 56 
and in crystallographic properties 57 . 

Our theoretical method and results to be presented 
here can be embedded in a path integral approach to the 
partition function of the underlying model. A part of the 
results was already discussed in 30 ' 5 . In the present paper 
we add further details of the formalism and present new 
results on the frequency dependence of the action kernel 
and the role of an external magnetic field. 

Moreover, we give a detailed discussion of the gradient 
of the carrier ground state energy with respect to fluctu- 
ations in the orientations of the localized Mn moments. 
For the case of isotropic systems with simple parabolic- 
band carriers the collinear ferromagnetic state with all 
Mn spins in parallel is always stationary (i.e. it has a 
vanishing energy gradient), but in general not stable 30 . 
In the present work we show the possibility of compla- 
nar spin structures which are local minima on the en- 
ergy landscape of such systems. For the case of models 
with spin-orbit anisotropy in the carrier system we find 
that the collinear state is not even stationary and does 
therefore not represent an energy minimum. This finding 
resembles recent results by Zarand and Janko 40 . 

Another study related to the present one was reported 
on recently by Korzhavyi et al. 4S . These authors inves- 
tigate, by both experimental and theoretical means, the 
influence of As antisite defects on the ground state struc- 
ture of the Mn local magnetic moments. As a result, 
As antisite defects are found to favor the formation of 
disordered local moment configurations of Mn spins with 
antiparallcl orientation and therefore reduce the net mag- 
netization. However, differently from the results of the 
present investigations, the Mn spins in the disordered lo- 
cal moment configurations are restricted to have (at ran- 
dom) either spin up or down with respect to some given 
quantization axis, i.e. they are collinear to each other. 

This paper is organized as follows. In section II we 
review our general theory of magnetic fluctuations in the 
kinetic exchange models studied here, and add further 
technical details. In section III we report on numeri- 
cal results concerning the instability of the collinear fer- 
romagnetic state in the parabolic-band model. In sec- 
tion IV we discuss several results arising from the very 
general gradient expression for the ground state energy 
obtained in section II. In subsection IV A 1 we demon- 
strate the occurrence of complanar or helical energy min- 
ima in the parabolic-band model. An important further 
finding is discussed in IV A 2 where the collinear state is 
shown not even to be stationary with respect to magnetic 
fluctuations if anisotropy induced by spin-orbit coupling 
is present. In subsection IV B we report on numerical 
steepest-descent studies of the true energy minima be- 
low the collinear state, and on the influence of an exter- 
nal magnetic field. In section V we discuss the small- 



frequency dependence of the effective-action kernel de- 
rived in II. We close with conclusions in section VI. 



II. MAGNETIC FLUCTUATIONS: GENERAL 
THEORY 

A. Path integral approach 

We now describe a path integral approach to the par- 
tition function of fluctuations around a given magnetic 
state of the Mn spins. Our method is similar to the one 
used in Refs. 25 with the differences that we do not ap- 
proximate the Mn magnetic moments by a continuum but 
retain them as individual localized spins, and we study 
fluctuations around more general noncollinear magnetic 
states. The state we consider is a tensor product of spin- 
coherent states for each Mn location I with a polarization 
direction 0/ = (sin #/ cos <pj, sin i?j sin <pj, cost?/). We 
introduce local coordinates for each Mn site whose in- 
direction coincides with 17/. The spin operators Tj with 
respect to these local bases are related to the original 
operators Si via 

Sf \ ( cos (pi — sin ifi \ 
Sj \ — \ sm <Pi cosiyS/ 

SfJ \ o 1 1 

( cost?/ sintfj \ ( Tf\ 
• 10 1*. (3) 

y -sint?/ cost?/ J \ Tf J 

To parametrize fluctuations around the magnetic state 
given by the directions 17/ we use the usual Holstein- 
Primakoff representation of the spin operators T/, 

T+ = ^2S-bjbibi , T z =S- b+bi . (4) 

Formulating the partition function as a path integral and 
representing the Mn spin Holstein-Primakoff bosons in 
a coherent-state parameterization one arrives after inte- 
grating out the carrier degrees of freedom at the following 
effective action: 

rP 

S cS = / dr V [zid T zi] - tr [In (G M V + SG^ 1 )] . (5) 
Jo / 

z/(t) stands for the bosonic Holstein-Primakoff field pa- 
rameterizing the fluctuations of the Mn spin I around 
its mean direction 17/. The integration over the imagi- 
nary time t goes from zero to the inverse temperature /?, 
and the trace in the second contribution is over fermionic 
carrier degrees of freedom and imaginary time. The the 
fluctuation-free (i.e. mean-field) part G~^ F of the inte- 
gral kernel is given by 
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where \x is a chemical potential and s z and s 1 * 1 = s x ± is y 
are carrier spin operators. In the case of a simple 
parabolic band they are just proportional to Pauli ma- 
trices while, for instance, in the case of a six-band k ■ p 
Hamiltonian they have a more complex form 54,26 . 

The fluctuation part of the inverse Green's function in 
up to quadratic order in the Holstein-Primakoff variables 
reads 5G' 1 = SG^ 1 + SG^ 1 with 



5G^ = W J ^-Ri) 



2 dj_ 
2 



- l <f"^2S [ zi cos" — -zi sin" — | s 



■ 2<M + 
2 



+ e +i *'V2S (z! cos 2 ^-zi sin 2 ^ s" 



- V25sini? 7 (z/ + z/)s 2 



(7) 
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(8) 



To analyze magnetic fluctuations we expand the action 
(5) in the bosonic spin variables describing deviations 
from the prescribed directions Cli, 

trln^V+.SG- 1 ) =trln(GW) 

oo 1 

-E-tr(-GMF*G-T. (9) 

n=l 

We will evaluate the effective action in up to second oder 
in the fluctuations zi(t). To this end we introduce the 

Fourier transforms Zi(Cl n ) = L dr exp(iQ n T)zi(T) with 
the notation zj(fl n ) = zi(Cl n )), where fi„ = 2nn/j3, n 
integer, is a bosonic Matsubara frequency. 

In this subsection we shall consider a simple parabolic- 
band model where fluctuation-free part of the carrier 
Green's function reads in real-space representation 



GMF^a;? \a';uj n ) = ~E 



iuj n ~ T\a 



(10) 



Here w„ = (2n+ l)-7r//3, n integer, is a fermionic Matsub- 
ara frequency and tp aa (r) is the spin component a of the 
carrier wave function with label a and energy e a = t] a +H- 



These wave functions are the eigenstates of the single- 



particle mean-field Hamiltonian TLmf = G 



MF 



The lowest-order contribution in the effective action 
occurs at zero Matsubara frequency only and is linear in 
the Holstein-Primakoff variables: 



with 5/ = g\ + igj, and 

g) = V2S(e VI x e z ) ■ J d 3 r[j(f- R T ) 

((W)) ■ e Vl )e Vl + ((s(r)) ■ e z ) e,) x Oj 



(11) 



gj = V2Se 



{e lf>1 x J 



d 3 rJ(r- Ri)(s{f}) 



(12) 
(13) 



Here (s(f)) is the expectation value of the carrier spin 
density, e Vl = (cos ipi , sin ipj , 0) and e z = (0,0, 1). 

(2) 

The contribution S^ c bilinear in the Holstein- 
Primakoff variables is quite complex for general Mn spin 
orientations 0/ . In particular, numerical evaluations of 
this quantity are extraordinarily tedious and computa- 
tionally expensive. We therefore shall concentrate on the 
collinear state where all Mn spins point, say, along the 
z-direction. Then the second-order contribution to the 
effective action reads 



41 = \ E E zi(n n )D u \n n )zAn n ) 



n I, J 



where the fluctuation matrix Djj(Q n ) reads 

DjJ(n n ) = Lu(n n ) + Ku(ci n ) 



with 



K 



u 



Su(- 

-IE 



riF{r) 

itt n +7] a -T]p 1 J 



a,0 I 

Here uf is the Fermi function, and 



(14) 

(15) 

(16) 
(17) 



J d 3 rJ( 



r-Ri)iJ a Ar)^^(r). (18) 



All quantities referring to the carrier system are to be 
evaluated for the collinear orientation of Mn spins. 

The diagonal contributions to the action kernel sum- 
marized in Ljj stem from the kinetic term Y]j zid T zi 
in the integrand of (5) and from the lowest-order term 
^(GmfSG^ 1 ) in (9). The term given by Kjj arise from 
the bubble contribution — (^{Gmf^G^Gmf^Gi 1 )) /2. 
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B. The static limit: Perturbation theory 

The limit of zero Matsubara frequency fi„ corresponds 
to static perturbations of the mean-field carrier ground 
state. It is instructive to verify this explicitly by elemen- 
tary perturbation theory. The unperturbed mean field 
Hamiltonian is again TLmf — G MF ~ d T + fi which is 
subject to a perturbation TLfi uc [{zi, Zi}] with 

TLfi uc [{z I ,z I }} = SG^[{z I ,z I }]+SG^ 1 [{z I ,z I }}. (19) 

Here 5G7/ 1 and SG^ 1 have the same form as in (7), (8) 



with zi = (T x + iT v )/\/2S parameterizing the pertur- 
bations with respect to the prescribed general directions 
Vtj in TLmf- 

We now consider the contributions from TLfi uc to the 
energy of the carrier ground state. In linear order in zj 
one has the contribution from the expectation value of 
SG^ 1 which can be written as 



(20) 



with gi = gj+igj given by (12), (13). This expression co- 
incides with (ll) 58 . The expectation values of the carrier 
spin density (s(r)) entering (20) via the complex coeffi- 
cients gi have to be computed within the ground state of 
the unperturbed Hamiltonian TLmf- This operator de- 
pends on the Mn spin directions fti and on the kinetic 
term TLkin for the carriers, which has not been specified 
so far within this perturbational approach to the ground 
state energy. Therefore the expression (20) holds for- 
mally for any band representation of the carrier system. 
In particular (20) is valid for the simple parabolic-band 
model as well as for more sophisticated k-p Hamiltonians. 
We will come back to this important fact in section III A. 

There are two contribution bilinear in zi . The first one 
is just the expectation value of SG^ 1 within the carrier 
ground state, 



4 2) = E *i l u*j 



i, j 



with 



-t„f 



d 3 rJ(r- fl/)(s(r))Q/ 



(21) 



(22) 



which obviously coincides with (16) for fi n = and fi/ = 
e z for all /. The other bilinear contribution is the second- 
order term arising from SG^ 1 which reads 



a, 13 



n F (Va) - n F (r)p) 



+ (cos ??/9?{z/} - i^{ Zl }) F«T,0i e -<vi 
+ (cos dM{ Zl } + &{ ZI }) F a ^e +ivi 



■ (23) 



It is straightforward to see that for in the collinear case, 
dj = e z for all /, this expression takes the form 



E k ] = E ^K U Zj 



(24) 



i, J 



with K u = K u (n n = 0) given by Eq. (17) 59 . 

We note that the perturbational approach described 
here is crucially different from the RKKY approximation 
often used in the theory of spin glasses 60 . There the cou- 
pling of the carrier spin density to the local moments 
is treated as a perturbation to the non-polarized carrier 
Fermi sea. This is justified provided that the free-carrier 
Fermi energy is large compared to the energy scale of 
the exchange coupling which is conveniently measured in 
terms of the carrier mean-field splitting A = JpdSNMm 
where Nmu is the density of Mn spins. However, for typ- 
ical parameters of (III,Mn)V systems the Fermi energy 
and A are fairly of the same order of magnitude with the 
former quantity being often even smaller. Therefore the 
free-carrier ground state is not a good starting point for 
perturbation theory. To explore the true ground state of 
the system described by the Hamiltonian (1) for param- 
eters realistic for (III,Mn)V semiconductors one should 
rather use the carrier ground state in the presence of a 
fully polarized Mn spin system as the starting point, as 
done here. 



III. INSTABILITY OF THE COLLINEAR 
FERROMAGNETIC STATE IN THE ISOTROPIC 
PARABOLIC-BAND MODEL 

A. General discussion of the energy gradient 

As already stressed in section II B, the expression (11) 
obtained for the ground state energy in lowest oder in the 
Holstein-Primakoff variables is very general, i.e. it is for- 
mally the same for any representation of the carrier sys- 
tem (e.g. parabolic band, six-band k -p Hamiltonian,...). 

Let us now examine further the coefficients given in 
Eqs. (12), (13) that enter the energy gradient (11). The 
imaginary part g\ is nonzero if and only if the vector 



mi 



d 3 rJ(r- Ri)(s(r)) 



(25) 



does not lie in the plane spanned by e Vl = 
(cos ip i , sin ipi , 0) and e z — (0,0,1). Similarly, the real 
part g\ is nonzero if and only if the projection of fhi 
onto the plane spanned by e Vl , e z is not collinear with 
the direction flj = (sin??/ cosy/, sin-di sinfj, cos??/) of 
the Mn spin at site /. In summary, for a given orientation 



4 



of the Mn local moments to be stationary in energy the 
local directions fli must be parallel or antiparallel with 
rfij at each Mn site I. An example for such a situation 
is the collinear ferromagnetic state with all Mn spins in 
parallel in a simple parabolic-band model. Here the spins 
of all carrier eigenstates are polarized along the common 
axis of the Mn spins, and therefore mj is always collinear 
with this direction at each site I. Thus, this state is al- 
ways stationary, but, as we shall see below, in general 
not stable. Moreover, the collinear ferromagnetic state 
turns out to be not even stationary if anisotropy induced 
by spin-orbit interactions is present, as it is the case for 
valence band holes in III-V semiconductors. 



B. Numerical results 

We now study the stability of the collinear ferromag- 
netic state in a simple parabolic-band model where the 
carriers are characterized by just one effective mass m*. 
As seen before, the energy gradient vanishes for such a 
state, and the quadratic zero-frequency (n = 0) contri- 
bution to the effective action (14) describes the energy of 
static fluctuations around the collinear state. Thus, for 
this state to be stable, the matrix Djj (0) must have non- 
negative eigenvalues only, while the occurrence of nega- 
tive eigenvalues of this matrix indicates that the perfectly 
collinear state is not the ground state. We note that 
for any arrangement of the Mn positions Ri the matrix 
Djj (0) contains a zero eigenvalue corresponding to a uni- 
form rotation of all spins. If Djj(0) is non-negative its 
eigenvalue distribution can be interpreted as a density of 
states (DOS) for magnetic excitations above the collinear 
state. 

We have evaluated the spectrum of Djj (0) in systems 
given by a simulation cube with periodic boundary con- 
ditions averaging over different realizations of the Mn 
positions. The single-particle wave functions il> aa (r) are 
computed in a plane- wave basis taking into account wave 
vectors q with length up to an appropriate cutoff q c . The 
same truncated plane- wave basis is used to compute the 
quantities (18) entering (17). Note that, for fluctuations 
around the collinear ferromagnetic state, Djjiiuj) is al- 
ways real and symmetric for real w since all carrier wave 
functions have for a given spin projection a a coordinate- 
independent phase (and can therefore also chosen to be 
real). This follows from the fact that the single-particle 
Hamiltonian describes for each spin projection just the 
problem of a spinless particle in a potential landscape 
provided by the Mn ions. Since Djj(ito) is real and sym- 
metric the components of each of its eigenvectors have all 
the same phase (and can be chosen to be real) . Physically 
this corresponds to the invariance of the system under 
rotations around the magnetization axis of the collinear 
state. 

The two upper panels of Fig. 1 show results for typi- 
cal system parameters for two different values of q c . The 



comparison of both panels shows that the effects of the 
wave vector cutoff on the low-lying excitations have al- 
ready saturated for the smaller q c . The by far largest 
contributions to the DOS lie at positive energies, while 
a small amount of sightly negative eigenvalues of Djj (0) 
indicate an instability of the perfectly collinear state. 

In the calculations discussed so far the Mn positions 
were chosen completely at random with uniform distri- 
bution, while in a real (III,Mn)V semiconductor the Mn 
ions are supposed to be located on the cation sites form- 
ing an fee lattice. In the bottom panel of Fig. 1 we show 
data for the same system parameters as in the top panel 
but with the Mn positions chosen from an appropriate 
fee lattice such that about 5 percent of all sites are oc- 
cupied. Both plots are practically identical indicating 
that our observations do not depend on this detail of the 
modeling. 

The shape of the eigenvalue distribution of the fluc- 
tuation matrix Djj(0) is quite sensitive in detail to the 
Mn density Nmu, the carrier density p, and the Hamil- 
tonian parameters m*, J p d, &o- In our numerics, we have 
extensively investigated this high-dimensional parameter 
space in regions realistic for (III,Mn)V semiconductors. 
The general finding is that the ground state of the sys- 
tem is generically noncollinear. To our observation there 
occur always negative eigenvalues of Djj(0) indicating 
the instability of the collinear state, provided that large 
enough wave vector cutoffs and system sizes are consid- 
ered. 

To analyze further the nature of this instability we 
consider the participation ratio 



p(E) 



NVJ2\ME)\ 4 



(26) 



where aj(E) is the 7-th component of the (normalized) 
eigenvector of Djj(0) with eigenvalue E, and the sum- 
mation goes over all NV Mn sites. This quantity is an 
estimate for the fraction of components of a(E) being 
substantially nonzero. For instance, if a vector contains 
exactly a fraction of p nonzero components of equal mod- 
ulus and all others being zero, its participation ratio is p. 
The largest participation ratio of unity is achieved for the 
zero-energy uniform rotation mode where all components 
of the corresponding eigenvector are equal. 

Fig. 2 shows the disorder-averaged participation ratio 
for the same situation as in the top panel of Fig. 1 . The 
negative-energy modes have clearly higher participation 
ratio than the eigenvectors at positive energy. This shows 
that the instability of collinear state is due to long-ranged 
dynamics involving a large fraction of the spins present 
in the system. Qualitatively the same observations arc 
made for other values of system parameters. 
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IV. NONCOLLINEAR GROUND STATES AND 
THE INFLUENCE OF A MAGNETIC FIELD 

In the previous section we have investigated the stabil- 
ity of the collinear ferromagnetic state in the parabolic- 
band kinetic-exchange model. This state is always sta- 
tionary (i.e. has a vanishing energy gradient) but not 
necessarily stable. In the present section we extend our 
ground state studies using the energy gradient expression 
(11). We will consider the case of simple parabolic bands 
as well as more sophisticated kinetic carrier Hamiltonians 
incorporating spin-orbit anisotropy. 

A. Stationary states 

1. Helical states in the parabolic-band model 

We now discuss a particular class of metastable states 
in the parabolic-band model. 

The energy gradient expression (11) can be employed 
in a numerical steepest-descent procedure to search for 
true energy minima. To this end, one starts with the 
collinear state with all spins pointing into the z-direction 
and steps down in energy by performing sufficiently small 
rotations of the Mn spins according to some negative 
eigenvector of Djj(0) (or some linear combination of 
them). Now consider the case that only one eigenvector 
is involved with all components having a certain common 
phase. The resulting orientations of Mn spins will all lie 
in the same plane spanned by the z-axis and a direction 
in the xy-plane that is determined by the above phase 
factor. Since the parabolic-band Hamiltonian is invari- 
ant under spin-flips with respect to this given plane, the 
local spin density (s(r)) and consequently the vectors rhj 
will also lie in this plane 61 . As seen in the previous sub- 
section, for such a situation the imaginary part gj of the 
coefficients vanishes identically for all I. When now ap- 
plying the energy gradient expression in a steepest-decent 
procedure the real part g\ of the gradient coefficients will 
just move the directions Qi within this plane, while the 
imaginary parts g\ remain strictly zero. Eventually this 
procedure will end up in a stationary state with all Mn 
spins lying in a plane prescribed by the initial departure 
from the collinear state. 

These conclusions are confirmed by explicit numerics 
where such complanar (or helical) stationary states are 
indeed observed. This states can be seen as bona fide 
local minima on the energy landscape. 

However, as explained in more detail in the follow- 
ing subsection, such type of energy minima occur only 
in isotropic models but not in systems with spin-orbit 
anisotropy. 



2. Nonstationarity of the collinear state in the presence of 
spin-orbit anisotropy 

We now examine the energy gradient (11) with re- 
spect to a six-band k ■ p Hamiltonian involving spin-orbit 
anisotropy 54,26,55 . In this case the rotational invariance 
in spin space is broken down to the cubic symmetry of 
the underlying GaAs crystal. 

Let us first consider the collinear ferromagnetic state. 
In this case we find numerically that the vectors fhi 
are not parallel (or antiparallel) with any given common 
orientation of the Mn spins. This finding includes also 
the crystallographic symmetry axes (1,0,0), (1,1,0), and 
(1,1,1) or their equivalents. 

There are two ways to restore the collinearity between 
the vectors rhi and a common orientation of the Mn 
spins: 

(i) Virtual crystal or continuum approximation: If the 
Mn spins are approximated as an continuum and point all 
along one of the crystallographic symmetry axes (1,0,0), 
(1,1,0), and (1,1,1) or their equivalents, the rhi are al- 
ways antiparallel to this direction. This finding explicitly 
confirms an assumption which earlier spin- wave calcula- 
tions in the six-band continuum model were based on 27 . 
Formally the continuum limit is reached by putting the 
parameter ao in the coupling function (2) to large values. 
In this limit the disorder with respect to the Mn positions 
is removed. 

(ii) Zero spin- orbit coupling: For vanishing spin-orbit 
coupling the vectors m; are antiparallel to any given com- 
mon direction of the Mn spins. 

Thus, in the presence of disorder with respect to the 
Mn positions and a finite anisotropy induced by spin- 
orbit interaction, the collinear ferromagnetic is never sta- 
tionary and therefore not the ground state. This state- 
ment includes the case that the common direction of the 
Mn spins is along some crystallographic symmetry axis. 

This important result resembles recent findings by 
Zarand and Janko 40 obtained within the RKKY approx- 
imation, where the ground state of the system was also 
found to be noncollinear (or, as termed there, frustrated). 

Finally we briefly remark on the case that the direc- 
tions of all Mn spins lie in some common plane. Accord- 
ing to our numerics, and as to be expected from the above 
results, the local spin density (s(r)) and consequently the 
vectors m/ do not lie in plane of the Mn spins. Therefore 
truly complanar states being local energy minima do not 
exist if spin-orbit anisotropy is present 

B. Steepest-descent results and the role of an 
external field in the parabolic-band model 

We have employed the energy gradient expression (11) 
in a numerical steepest-descent procedure outlined in 
subsection IV A 1 to search for stationary states in the 
parabolic-band model. These states can be considered as 



G 



bona fide energy minima. Our results are as follows: In 
cases where the energy minimum found by this method 
is close to the collinear state (with a magnetization of 
about 90 percent of the maximum value or more), this 
minimum appears to be unique (for a given disorder real- 
ization) and can therefore be considered as the true abso- 
lute ground state of the system. However, in situations 
where the magnetization is reduced more substantially 
(by, say, about 20 percent or more) the energy minima 
found are not unique anymore but depend on technical 
details of the numerical procedure. In such cases the sys- 
tem is essentially spin-glass like with a complicated en- 
ergy landscape. This situation occurs typically at larger 
density ratios p/Nmu of carriers and Mn spins. For the 
system shown in Fig. 1 for instance the magnetization 
values in the energy minima typically found are of about 
30 to 40 percent of the collinear state. 

Finally we consider the role of a magnetic field coupled 
to the Mn spins, 



Hb = A 



E^ 



(27) 



where A — gfisSB is the vector of Zeeman couplings. 
Here the magnetic field is coupled only to the Mn spins 
and not to the band carriers, since their mean-field spin 
splitting A = J p dSN Mn is large compared to the Zeeman 
couplings considered here. Therefore this contribution to 
the effective carrier spin splitting is negligible, and one 
should expect the collinear state to be stabilized by an 
external Zeeman field which is equal in magnitude to the 
smallest (negative) eigenvalue of the fluctuation matrix 
Djj (0). For instance, from the data shown in Fig 1 this 
Zeeman field would be about 0.5meV. 

A magnetic field in the z-direction leads to an ad- 
ditional contribution to the gradient components g\ of 
A z -\/2/5sini?/, which we shall consider in the following. 
We employ the following numerical procedure. Starting 
from the collinear state with all Mn spins pointing in the 
z-direction we obtain a noncollinear energy minimum by 
the steepest-descent method described above. Then we 
add a magnetic field in the z-direction and an repeat the 
procedure to reach a new energy minimum. This state is 
then the new starting point for a steepest-descent walk 
with an increased magnetic field. By iterating this proce- 
dure we obtain a zero-temperature magnetization curve 
as a function of the external field. 

Our numerical results are as follows: In cases where 
the zero-field ground state magnetization is substantially 
reduced from its maximum value for the collinear state 
(i.e. in the "glassy" regime), the magnetization at fi- 
nite field depends quantitatively on the step width that 
is used in increasing the field. Fig. 3 shows magnetiza- 
tion data for the same system parameters as in the data 
before. As seen from there the Zeeman field needed to 
realign all spin along the z-direction is somewhat larger 
than the expected value of A z — 0.5meV, which is a 
typical hysteresis effect. 



V. THE SMALL-FREQUENCY DEPENDENCE 
OF THE ACTION KERNEL AND ADIABATIC 
DYNAMICS 

In section III we studied the stability of the collinear 
ferromagnetic state within a parabolic band model in 
terms of the static component of the effective action ker- 
nel Djj. Here we extend these analysis to the regime of 
small but finite Matsubara frequencies Q n . 

Expanding Ku(fl n ) around £1„ = one has 



Ku(Cl n ) = Ku(0) + ifl n K' u + 



with 



K' 



ij 



a, /3 



n F (y a ) - n F (r]p) 



(28) 



(29) 



The matrix K' fulfills an important sum rule which 
we derive now. As mentioned before, for the collinear 
state with all Mn spins pointing in the z direction all 
carrier eigenstates have cither spin up or spin down, 
and one ends up with two separated problems for spin- 
less fermions in a potential landscape. The stationary 
Schrodinger equations for spin down and spin up parti- 
cles read 



2m 



i J 

■C^T^ + f E / d 3 rJ(r-Rj)^(r} 



(30) 



= e^ T (f). (31) 

Taking all wave functions to be real and combining these 
equations one finds 



a|,/3T 



(32) 



where (-|-) denotes a scalar product between spinless wave 
functions. Since both the spin up and the spin down 
carrier wave functions fully span the Hilbert space of a 
single spinless particle we have 



S 2 



1 



= E 



{W - Va)' 



IJ 



and similarly 



E 



ol (Vl3-V»f j j 1 J 

From these sum rules one derives 62 



(33) 



(34) 
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where n-\, ri| are the numbers of up and down spin car- 
riers, respectively. 

We now consider the limit of large inverse temperature 
(3 where the bosonic Matsubara frequencies fi n = 2nn/(3 
can be treated as a continuous variable ft. We are inter- 
ested in the adiabatic limit of the spin dynamics which 
is described by the low-frequency limit of the effective 
action (14). In the expansion of the Fourier transformed 
Holstein-Primakoff variables, 

zi(fi) = zi(O) + fi(J0 +••• (36) 

the zero- frequency component 2/(0) corresponds to 
Holstein-Primakoff variables constant in imaginary time 
r. Here we consider the uniform case Zi(t) = 
\fSrj exp(«x) for all /. This describes a uniform rotation 
of all Mn spins around the axis (sinx, — cosx, 0) by an 
angle parametrized by n. As seen in section III the ma- 
trix Djj(0) annihilates the vector of the zero-frequency 
components zi(Q = 0) = fly/Srj cxp(ix)- Therefore, us- 
ing the sum rule (35), the low- frequency expansion of the 
effective action (14) reads 

4uc = m)v (-SN Mn V + \ (n T - +■■■ (37) 

where NmuV is the number of Mn ions in the system. 
With an analytical continuation to real time, the above 
expression linear in O is the geometric phase generated 
by an adiabatic uniform rotation of all Mn spins (point- 
ing initially along the ^-direction) around an axis in the 
xy-plane by an angle /3J7r/. The second term in the paren- 
theses stems from the carriers which adiabatically follow 
the Mn spins. Note that this contribution comes with a 
different sign. This is due to the fact that we are deal- 
ing with an effective action for the Mn spins only where 
the carriers have been integrated out, and therefore only 
the Mn spins are 'actively' rotated. In our formalism the 
initial polarization axis of the Mn spins defines the quan- 
tization axis for the carrier spins, and rotating this axis 
is just a 'passive' rotation of the carrier spin coordinate 
system. Therefore the geometric phase stemming from 
the carriers has a different sign. 

We now turn to the case of general (i.e. non-uniform) 
adiabatic rotations of the Mn spins parametrized by other 
eigenvectors of Dyj(0) than the one discussed above. 
Then the geometric phase of the localized spins still stems 
from the frequency-dependent part of the diagonal ma- 
trix Ljj(O) (cf. Eq. (16)) while the largest contribu- 
tion to the carrier phase will arise from Djj(0) itself 
and the derivative K'jj. We therefore expect K\j to be 
dominantly diagonal when expressed in the eigenbasis of 
Djj(0). In fact, this expectation is confirmed by numer- 
ical evaluations of this quantity. As a measure for the 
'diagonality' of a given matrix A we consider 



D{A) ^^§-) 1/2 ■ (38) 

This quantity is unity if A is diagonal, and of order the 
inverse of the square root of the dimension of A if all of 
its elements are of the same order of magnitude. 

We have evaluated the 'diagonality' D of K\ } ex- 
pressed in the eigenbasis of Djj(0) for various sets of 
system parameters and averaged this quantity over many 
disorder realizations with respect to the Mn positions in 
space. We consistently find D substantially larger than 
the average value of a general matrix. For instance, for 
the parameters used in Figs. 1,2 we find D = 0.560 (av- 
erage over all disorder realizations) with a fluctuation 

AD = v D 2 — D 2 = 0.045. This is by an order of mag- 
nitude larger than the value of a general matrix of this 
dimension which is D = 1/V400 = 0.05. 

VI. CONCLUSIONS 

We have studied the ground state properties of kinetic- 
exchange models for carrier-induced ferromagnetism in 
(III,Mn)V semiconductors with randomly distributed Mn 
ions. Our method is embedded in a path integral spin- 
wave type formalism leading to an effective action for 
the Mn spins with full Matsubara frequency dependence. 
The zero-frequency contribution to this action is equiv- 
alent to static perturbation theory and characterizes the 
stability of a given spin configuration, while the compo- 
nent linear in frequency can be interpreted as the joint 
Berry phase of the Mn and carrier system. 

Our perturbational approach to the ground state of 
the system studied here differs from the RKKY approx- 
imation in so far as we not do do perturbation the- 
ory around the free-carrier ground state but around the 
carrier ground state in the presence of a fully aligned 
collincar Mn spin system. This is appropriate since the 
mean-field spin splitting is for realistic system parame- 
ters not small compared to the Fermi energy. Therefore 
it cannot be regarded as a small perturbation to the free- 
carrier ground state, as done in the RKKY approxima- 
tion. This approximation actually works well for metallic 
spin glass systems where the coupling to local moments 
is indeed a small perturbation 60 . 

For parabolic-band carriers the collinear ferromagnetic 
state with all Mn spins in parallel is always stationary 
but generically unstable. This instability can be charac- 
terized in terms of inverse participation ratios and is due 
to long-ranged nonlocal spin fluctuations. We also have 
presented results for the ground state magnetization as 
a function of an external field. 

For carrier dispersions involving anisotropy induced by 
spin-orbit coupling the collincar state is not even station- 
ary and therefore also not the ground state. This inter- 
play between the anisotropy in the carrier system and 
the disorder in the Mn positions reflects recent findings 
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by Zarand and Janko 40 obtained within the RKKY ap- 
proximation. The stationarity of the collinear state is 
restored in the continuum or virtual crystal approxima- 
tion where disorder is neglected. 
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FIG. 1. The disorder-averaged density of states of magnetic 
excitations for a simulation cube of volume V = L = 400nm 
with a Mn density of Nmu = 1.0nm~ 3 and a density of 
p = 0.15nm~ 3 of carriers having a band mass of half the 
bare electron mass. The strength of the exchange interaction 
between ions and carriers is J pc i = 0.05eVnm~ 3 with a spatial 
range of ao = 0.40nm The two upper panels show data for 
different wave vector cutoff q c with the Mn positions chosen 
completely at random. The lowest panel contains data for 
the same situation as the top one but with the Mn positions 
chosen from an fee lattice. The peaks at zero energy are due 
to the uniform rotation mode which strictly occurs in any 
disorder realization. 
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value at zero energy is enhanced due to the contribution of 
the uniform rotation mode in that sample interval. 
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FIG. 3. The ^-component of the Mn spin polarization 
£l z = y^ r flj I (NmuV) as a function of an external magnetic 
field for the same parameters as in the top panel of Fig. 1. 
The Zeeman coupling is increased in units O.lmeV. The thin 
lines are results for three individual disorder realizations while 
the thick graph is the disorder average over 11 realizations. 
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